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ABSTRACT: Vacuum polarisation in QED in a background gravi- 
tational field induces interactions which effectively violate the strong 
equivalence principle and affect the propagation of light. In the low fre- 
quency limit, Drummond and Hathrell have shown that this mechanism 
leads to super luminal photon velocities. To confront this phenomenon 
with causality, however, it is necessary to extend the calculation of the 
phase velocity f p h(w) to high frequencies, since it is f p h(oo) which de- 
termines the characteristics of the effective wave equation and thus the 
causal structure. In this paper, we use a recently constructed expres- 
sion, valid to all orders in a derivative expansion, for the effective action 
of QED in curved spacetime to determine the frequency dependence 
of the phase velocity and investigate whether superluminal velocities 
indeed persist in the high frequency limit. 



SWAT/333 
March 2002 



1. Introduction 



It has been known since the original work of Drummond and Hathrell 0] that quantum 
effects have important consequences for the propagation of light in curved spacetime. 
In the classical theory of electrodynamics coupled to general relativity, light propagates 
simply along null geodesies. In quantum electrodynamics, however, vacuum polarisation 
changes the picture and the background gravitational field becomes a dispersive medium 
for the propagation of photons. In itself, this is perhaps not surprising. The one-loop 
vacuum polarisation contribution to the photon propagator introduces a non-trivial length 
scale A c (the inverse electron mass) and it is natural that photon propagation will be 
significantly affected when the typical curvature scale L is comparable to A c . However, 
the remarkable result found by Drummond and Hathrell is that in many cases the effect of 
vacuum polarisation is to induce a change in the velocity of light to 'superluminal' speeds, 
i.e. v > c, where c is the fundamental constant. 

Since this original discovery, many special cases have been studied in detail 
including propagation in black hole spacetimes described by the Schwarzschild, Reissner- 
Nordstrom or Kerr metrics, the FRW metric of big bang cosmology, and gravitational 
wave backgrounds, in particular the Bondi-Sachs metric describing asymptotic radiation 
from an isolated sourceB. The phenomenon of superluminal propagation has been observed 
in all these examples of gravitational fields and a number of general features have been 
identified. A notable result is the 'horizon' [H] or 'touching' [15] theorem, which shows 



that even in the presence of superluminal velocities, the effective black hole event horizon 
defined by physical photon propagation coincides precisely with the geometric horizon. 
Another important observation is that the speed of light increases rapidly (as 1/t 2 in 
the radiation dominated era) in the early stages of a FRW cosmology [IJ, with potential 
implications for the horizon problem and related issues. 

All this work has been based on the initial Drummond-Hathrell analysis, in which they 
showed that the effect of one-loop vacuum polarisation is to induce the following effective 
action, generalising the free-field Maxwell theory: 



r = / dx\f^g 



+ -L (aRF^F^ + bR^F^F'x + cR^ Xp F^ F Xp + dD^ x D„F u X ^J 



(1.1) 



1 Other examples and a selection of related work may be found in refs.§,|,0,|,|,|O) and 
possible implications for time machines in [11, 12, 13 1 . 
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Here, a,b,c,d are constants of 0(a) and m is the electron mass. The notable feature is the 
direct coupling of the electromagnetic field to the curvature. This is an effective violation 
of the strong equivalence principle, which states that dynamical laws should be the same 
in the local inertial frames at each point in spacetime. More precisely, this requires the 
coupling of electromagnetism to gravity to be through the connections only, independent 
of curvature. Eq.(TTT) shows that while this principle may be consistently imposed at the 



classical level, it is necessarily violated in quantum electrodynamics. 

Standard geometric optics methods applied to this effective action results in a modified 
light cone |I|,^4[: 

k 2 - %(2b + Ac)T^k x + ^-C^ Xp k»k x a v a p = (1.2) 



where k p is the wave vector and a M is the polarisation vector. In the second term, we have 
replaced the Ricci tensor by the energy- momentum tensor using the Einstein equations. 
This emphasises that this contribution is related to the presence of matter; indeed this 
takes the same form as the subluminal corrections to the speed of light in other scenarios, 
such as background magnetic fields, finite temperature, etc. [ |16l , |T7|| . The uniquely gravi- 
tational term involving the Weyl tensor depends explicitly on and gives a polarisation 
dependence of the speed of light (gravitational birefringence). Moreover, it is readily seen 
that this term changes sign for the two physical, transverse polarisations so that for vacuum 
(Ricci flat) spacetimes, if one photon polarisation is subluminal, the other is necessarily 
superluminal. This property, and many others, are most clearly seen by rewriting the light 
cone condition (|1.2|) in Newman- Penrose formalism |14| . 



The effective action ( p. . 1|) is, however, only the first term in a derivative expansion, 
with higher order terms in 0( — ) omitted. The corresponding light-cone condition (1.2) is 



therefore valid a priori only in the low frequency approximation. The modified light 
velocity derived from it is the phase velocity i> p h = ^ (where k^ = (u,k) in a local inertial 
frame) at uj ~ 0. In order to discuss the obvious issues concerning causality, however, the 
relevant 'speed of light' is not i> p h at u ~ but t> p h in the high frequency limit u — > oo. 
(We discuss this point carefully in section 2.) In order to address causality, therefore, we 
first need to establish the effective light cone condition for high frequencies, which itself 
involves finding the 'high frequency limit' of the effective action (|1 . 1|) . 

We have recently evaluated the effective action for QED in curved spacetime to all 

orders in 0( — ), keeping terms of the form RFF as in (|1 . 1| ) , that is, the terms relevant to 

r a 2 

photon propagation to lowest order in O(-^j), i.e. 0(j%). The derivation and full details 



of this result are presented in ref. |18[ . Here, we generalise the geometric optics derivation 
of photon propagation using this new effective action and discuss what we can learn about 
the critical high frequency behaviour of the phase velocity. The question we wish to answer 
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is whether the phenomenon of superluminal photon propagation is a curiosity of the low 
frequency approximation or whether it persists at high frequency, forcing us to confront 
the serious implications for causality associated with faster than light motion. In the 
end, our results are intriguing but not as yet conclusive. It appears that further field- 
theoretic developments may be needed to give a final resolution of the nature of dispersion 
in gravitational fields. 

The paper is presented as follows. In order to clarify exactly what we mean by the 
"speed of light" and why it is u p h(oo) which controls the causal behaviour of the theory, 
we review various definitions and theorems concerning the propagation of light in section 
2. In section 3, we review the fundamentals of geometric optics and its application to 
the Drummond-Hathrell action, clarifying some subtle points arising from earlier work. 
Our new result for the QED effective action in curved spacetime valid to all orders in the 
derivative expansion is presented in section 4, including some technical formulae for form 
factors. The implications for photon propagation in the high frequency limit are described 
in sections 5 and 6, leading to the apparent prediction that high frequency superluminal 
velocities are possible in certain spacetimes. This preliminary conclusion is challenged in 
section 7, where we compare our gravitational analysis with the closely related problem 
of photon propagation in a background magnetic field. Section 8 summarises our final 
conclusions. 



2. The "Speed of Light" 

Our fundamental interest is in whether the fact that the phase velocity at low fre- 
quencies can be superluminal (and therefore imply motion backwards in time in a class 
of local inertial frames) is in contradiction with established notions of causality. It has 
been discussed elsewhere that while in special relativityi superluminal motion necessarily 
implies a causal paradox, in general relativity this is not necessarily so [PUl4fl- The key 
question is whether a spacetime which is stably causal (see ref. [[23], Proposition 6.4.9 for a 
precise definition) with respect to the original metric remains stably causal with respect 
to the effective metric defined by the modified light cones. Essentially, this means that the 
spacetime should still admit a foliation into a set of hypersurfaces which are spacelike ac- 
cording to the effective light cones. This is an interesting global question, which is beyond 
the scope of this paper. It seems entirely possible, however, that the Drummond-Hathrell 



2 To be precise we mean here: "in Minkowski spacetime with no boundaries" . An analysis of 
causality related to the phenomenon of superluminal propagation between Casimir plates []l9| , p0| 

has recently been given in ref. (see also nnm). 
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modifications to the light cones should not destroy stable causality. 



Before addressing such issues, however, we need to be clear what exactly we mean by 
the "speed of light" which determines the light cones to be used in determining the causal 
structure of spacetime. In this section, we therefore review briefly some basic definitions 
and results from classical optics in order to motivate our subsequent analysis. 

A particularly illuminating discussion of wave propagation in a simple dispersive 
medium is given in the classic work by Brillouin |P4"| . This considers propagation of a 
sharp-fronted pulse of waves in a medium with a single absorption band, with refractive 
index n(u): 



n 2 (u) = 1 



a 



(2.1) 



u 2 — Uq + 2iujp 

where a,p are constants and c„ is the characteristic frequency of the medium. FivJ 
distinct velocities are identified: the phase velocity v p h = Cjj^ = K^jy, group velocity 



dui 



f gp = ^pry, signal velocity t> S i g , energy-transfer velocity v en and wavefront velocity v w {, with 
precise definitions related to the behaviour of contours and saddle points in the relevant 
Fourier integrals in the complex w-plane. Their frequency dependence is illustrated in 
Fig. 1. 
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Fig.l Sketch of the behaviour of the phase, group and signal velocities with frequency in the model 
described by the refractive index (2.1). The energy-transfer velocity (not shown) is always less than c and 
becomes small near ujq. The wavefront speed is identically equal to c. 



3 In fact, if we take into account the distinction discussed in section 3 between the phase 
velocity and the ray velocity v ray , and include the fundamental speed of light constant c from the 
Lorentz transformations, we arrive at 7 distinct definitions of "speed of light" . 
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As the pulse propagates, the first disturbances to arrive are very small amplitude waves, 
'frontrunners', which define the wavefront velocity t> w f- These are followed continuously 
by waves with amplitudes comparable to the initial pulse; the arrival of this part of the 
complete waveform is identified in [[24[ as the signal velocity v s { g . As can be seen from 
Fig. 1, it essentially coincides with the more familiar group velocity for frequencies far 
from ujq, but gives a much more intuitively reasonable sense of the propagation of a signal 
than the group velocity, whose behaviour in the vicinity of an absorption band is relatively 
eccentric. @ As the figure also makes clear, the phase velocity itself simply does not represent 
a 'speed of light' relevant for considerations of signal propagation or causality. 

The appropriate velocity to define light cones and causality is in fact the wavefront 
velocity t> w f. (Notice that in Fig. 1, v w f is a constant, equal to c, independent of the 
frequency or details of the absorption band.) This is determined by the boundary between 
the regions of zero and non-zero disturbance (more generally, a discontinuity in the first 
or higher derivative of the disturbance field) as the pulse propagates. Mathematically, 
this definition of wavefront is identified with the characteristics of the partial differential 



equation governing the wave propagation ||27|1 . Our problem is therefore to determine the 



velocity associated with the characteristics of the wave operator derived from the modified 
Maxwell equations of motion appropriate to the new effective action. 

A very complete and rigorous discussion of the wave equation in curved spacetime is 
given in the monograph by Friedlander |[28| , in which it is proved (Theorem 3.2.1) that 



the characteristics are simply the null hypersurfaces of the spacetime manifold, in other 
words that the wavefront always propagates with the fundamental speed c. However, 
this discussion assumes the standard form of the (gauge-fixed) Maxwell wave equation 



(cf. ref. ^8|, eq. (3.2.1)) and explicitly does not cover the modified wave equation derived 



from the action (|1.1|) , precisely because of the extra curvature couplings. 

Instead, the key result for our purposes, which allows a derivation of the wavefront 
velocity, is derived by Leontovich pi. In this paper, an elegant proof is presented for a 



very general set of PDEs that the wavefront velocity associated with the characteristics is 
identical to the u —* oo limit of the phase velocity, i.e. 

i> w f = lim — = lim v ph (u) (2.2) 

\rv\ 



This proof appears to be of sufficient generality to apply to our discussion of photon 
propagation using the modified effective action of section 4. 



4 Notice that it is the group velocity which is measured in quantum optics experiments which 
find light speeds of essentially zero |25|1 or many times c [pq|. A particularly clear description in 



terms of the effective refractive index is given in [25|. 
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The wavefront velocity in a gravitational background is therefore not given a priori 
by c. Taking vacuum polarisation into account, there is no simple non-dispersive medium 
corresponding to the vacuum of classical Maxwell theory in which the phase velocity repre- 
sents a true speed of propagation; in curved spacetime QED, even the vacuum is dispersive. 
In order to discuss causality, we therefore have to extend the original Drummond-Hathrell 
results for v p h(u ~ 0) to the high frequency limit v p h(uj — > oo), as already emphasised in 
ref. . This is why the effective action ( |1 . 1|) to lowest order in the derivative expansion is 
not sufficient and we require the all-orders effective action of section 4. 

A final twist emerges if we write the standard dispersion relation for the refractive 
index n{oS) in the limit uj — > oo: 



n( oo 



o r°° a, , 

n(0)-- / — $Sn(u) (2.3) 

7T Jo W 



For a standard dispersive medium, $sn(u) > 0, which implies that n(oo) < n(0), or equiv- 
alently t> p h(oo) > f p h(0). Evidently this is satisfied by Fig. 1. The key question though 
is whether the usual assumption of positivity of $sn(u) holds in the present situation of 
the QED vacuum in a gravitational field. If so, then (as already noted in ref.|I[]) the su- 
perluminal Drummond-Hathrell results for i> p h(0) would actually be lower bounds on the 
all-important wavefront velocity i> p h(oo). However, it is not so clear that positivity of 
3n(w) is reliable in the gravitational context. Indeed it has been explicitly criticised by 
Dolgov and Khriplovich in refs.|3(],|31|, who point out that since gravity is an inhomoge- 



neous medium in which beam focusing as well as diverging can happen (see next section), 
a growth in amplitude corresponding to $sn(u) < is possible. It therefore seems best to 
set the dispersion relation (|2.3j ) aside for the moment and concentrate instead on a direct 
attempt to determine w p h(oo). This is the goal of this paper. 



3. Low Frequency Photon Propagation 

The simplest way to deduce the light-cone condition implied by the effective action is 
to use geometric optics. In this section, we review the approach introduced in ref.[[l] (see 
also refs.p],[rS|) emphasising some points which will be important later. 



In geometric optics (see [[52J for a thorough discussion) the electromagnetic field is 
written as the product of a slowly-varying amplitude and a rapidly-varying phase, i.e. 

F liV = ft{f liV +ieh IJU , + ...)e>* (3.1) 

Here, e is a parameter introduced purely as a book-keeping device to keep track of the 
relative order of magnitude of terms. The field equations and Bianchi identities are solved 
order by order in e. 
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The wave vector is defined as the derivative of the phase, i.e. = <9 M $. The leading 

order term in the Bianchi identity D^F^ = is of 0{\) and constrains to have the 
form 

f/j,u = k^Av k v A^ (3-2) 

where A^ = Aa^. A represents the amplitude while (normalised so that a^a^ = —1) 
specifies the polarisation. For physical polarisations, k^a^ = 0. 

Conventional curved spacetime QED is based on the usual Maxwell action, so the 
equation of motion is simply 

D^F^ = (3.3) 

At leading order, 0(|), this becomes 

kj^ = (3.4) 

and since this implies 

fcV = (3.5) 

we immediately deduce that k 2 = 0, i.e. k^ is a null vector. Then, from its definition as a 
gradient, we see 

k^D^k u = k il D v k il = \D u k 2 = (3.6) 

Light rays (photon trajectories) are defined as the integral curves of i.e. the curves 
x M (s) where = These curves therefore satisfy 

o = fc , B ^ = _ + r^— — (3.7) 

which is the geodesic equation. So in the conventional theory, light rays are null geodesies. 
The subleading, O(l), term in the equation of motion gives 

WD^A V = - 1 -{D^)A v (3.8) 

which decomposes into 

k^DX = (3.9) 

and 

k^D^{\nA) = -^D^ (3.10) 

The first shows that the polarisation vector is parallel transported along the null geodesic 
rays while the second, whose r.h.s. is simply minus the optical scalar 6>, shows how the 
amplitude varies as the beam of rays focuses or diverges. 



We now apply the same methods to the modified effective action ( |1 . 1| ) ■ This gives rise 
to a new equation of motion which, under the approximations listed below, simplifies to: 



n - 

A 4 2 



2bR tlX D»F Xv + AcR u v Xo D^F xp 



(3.11) 



Here, we have neglected derivatives of the curvature tensor, which would be suppressed 
by powers of 0(X/L), where A is the photon wavelength and L is a typical curvature 
scale, and we have omitted the new contributions involving D^F^ : since this term is 
already 0(a) using the equations of motion, these contributions only affect the light cone 
condition at 0(a 2 ) and must be dropped for consistency Making the standard geometric 
optics assumptions described above, we then find the new light cone condition: 



1 



77r 



2bW x kJ Xv + 4cR^ Xp kJ Xp 







(3.12) 



Eq.( p.l2|) can now be rewritten as an equation for the polarisation vector a M , and re- 
expressing in terms of the Weyl tensor we find 



fcV - ( 2b + 4c ) R x (k»k x a v - k»k v a x ) - ^C» Xp k»k x a 
m z m 2 







(3.13) 



The solutions of this equation describe the propagation for a photon of wave vector /c M and 
polarisation a M . Contracting with a v , we find the effective light cone 



k" 



(2b + 4c) 



8c 



m 



Ra\k^k x + — ttC ' av\pk^k x a h a p 







(3.14) 



from which (|1.2|) follows immediately. 



Notice that in the discussion of the free Maxwell theory, we did not need to distinguish 
between the photon momentum i.e. the tangent vector to the light rays, and the wave 
vector hp since they were simply related by raising the index using the spacetime metric, 
pH _ g^ v ]^ v- in the modified theory, there is an important distinction. The wave vector, 
defined as the derivative of the phase, is a covariant vector or 1-form, whereas the photon 
momentum/ tangent vector to the rays is a true contravariant vector. The relation is 
non-trivial. In fact, when as in ( |3.14| ) we can write the light cone condition for the wave 
vector as the homogeneous form 

G^k^K = (3.15) 
we should define the corresponding 'momentum' as 



P 



(3.16) 
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and the light rays as curves x M (s) where ^j— = p^ . This definition of momentum satisfies 

G^tf = Q^k^K = (3.17) 

where G = Q~ x therefore defines a new effective metric which determines light cones 
mapped out by the geometric optics light rays. (Indices are always raised or lowered using 
the true metric g^v) The ray velocity v ray corresponding to the momentum p M , which is 
the velocity with which the equal-phase surfaces advance, is given by 

\E\ d k\ /oioN 

tW = -o=^f ( 3 - 18 ) 
along the ray, and is in general different from the phase velocity 

£.0 

v P h = (3.19) 



A nice example of this is given in ref. |13|], which analyses certain aspects of superlu- 



minal propagation in Minkowski spacetime with Casimir plates. The discrepancy between 
f p h and f ray (called cn g ht in jTSJ) in that example is due to a difference between the di- 
rection of propagation along the rays and the wave 3-vector. Otherwise, it follows directly 
from (|3.17| ) that v ray and v p h are identical. 

Recognising the distinction between p and k also clarifies a potentially confusing point 
in the important example of propagation in a FRW spacetime pj. Since the FRW metric 
is Weyl flat, the modified light cone condition (|1.2|) reads simply 

k 2 = ( T^k v (3.20) 

where £ = ^(2b + Ac) and the energy-momentum tensor is 

Tp V = (p + P)n^n v - Pg^ (3.21) 

with = (et)^ specifying the time direction in a comoving orthonormal frame, p is 
the energy density and P is the pressure, which in a radiation-dominated era are related 
by p — 3P = 0. The phase velocity is independent of polarisation and is found to be 
superluminal: 

k° 1 

v ph = — = l + -((p + P) (3.22) 

At first sight, this looks surprising given that k 2 > 0. However, if instead we consider the 
momentum along the rays, we find 

p 2 = g^v u = -Op + P)(k ) 2 (3.23) 
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and 

v ray = l 4 = l + U(p + P) (3.24) 
p z 

The effective metric G = Q~ l is 



n + C{p + p) 








\ 





-1 














-1 





V 








-1/ 



(3.25) 



In this case, therefore, we find equal, superluminal velocities t> p h = v ray and p 2 < is 
manifestly spacelike as required. 

In the radiation dominated era, where pit) = 3§^~ 2 , we have 

v ph = l + t~ 2 (3.26) 

which, as already observed in increases towards the early universe. Although this 
expression is only reliable in the perturbative regime where the correction term is small, 
it is intriguing that QED predicts a rise in the speed of light in the early universe. It 
is a matter of speculation whether this superluminal effect persists for high curvatures 
(L ~ A c ) and whether it could be important in the context of the horizon problem. 



4. The Effective Action for QED in a Gravitational Field 

The effective action presented in this section has recently been derived [[18| by adapting 
the more general background field calculations of []33|] , this latter paper being the culmina- 



tion of the theoretical development described in . The result is an effective action which 
incorporates the one-loop vacuum polarisation contributions to the photon propagator in 
an arbitrary, weak gravitational field. We therefore keep terms of type RFF, i.e. quadratic 
in the electromagnetic field but only of first order in the curvature. This neglects terms of 

A 2 

higher orders in 0{j% ). However, the new feature compared to the Drummond-Hathrell 



action ( p. . 1|) is that terms to all orders in derivatives are kept. This allows a discussion of 
the frequency dependence of the modifications to photon propagation. 

The effective action is given by 

T = r (0) +lndet^(x,a; / ) (4.1) 

where T(o) is the free Maxwell action and S(x, x') is the Green function of the Dirac 
operator in the background gravitational field, i.e. 



(ip-m)S(x,x') = -?=5(x,x') (4.2) 
V~9 
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In fact it is more convenient to introduce the scalar Green function G(x, x') defined 



by 



so that 



S(x, x') = (ip + m)G(x, x') 
1 



[D 2 + iea pv F^ -^R + m 2 ) G(x, x') 



-g 



-.5(x, x') 



(4.3) 
(4.4) 



Then we evaluate Y from the heat kernel, or proper time, representation 

r = r (0) - - / -e-™ s T?g{x,x'-s) 

1 Jo s 

where 



d 

Vg(x, x ; s) = i—Q(x, x ; s) 

OS 

with Q{x, x'; 0) = G(x, x'). Here, V is the differential operator in eq.(fOj) at m = 



(4.5) 



(4.6) 



The details of the derivation of the effective action are given in ref. [T^] and here we 
simply quote the result: 



dx\J —g 



-F^F^ + -^(D p F pX GlD v F» x 
4 m z V 

+ Gl RF^F"" + Gl R P „F» X F V X + Gl R pvXp F^F Xp 

+ -^-(gIrd^f pX d u f v x 

+ G~l R^D x F Xp D p F pv + Gl R pv D»F Xp D v F Xp 
+ G* 7 R^D^D v F Xp F Xp + Gl R pv D»D x F Xp F pv 

+ Gl R PiVXp D a F ap D x F^ u ) 

(4.7) 



In this formula, the G n [n > 1) are form factor functions of three operators: 



G n = G r 



n 2 n 2 d 2 

U {\) U (2) U (3) 



m z m z m z 



where the first entry (D 2 ^) acts on the first following term (the curvature), etc. Go is 
similarly defined as a single variable function. The G n are themselves expressed as proper 
time integrals: 



G n (x 1 ,x 2 ,x 3 ) 



1 a 
2n 



ds 



e ls (is) p g n (—isxi, —isx2, —isxs) 



(4.9) 
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where p = 1 for n = 0, . . . 3 and p = 2 for n = 4, . . . 9, and we have rescaled s by a factor 
of m? so as to be a dimensionless variable. 



A crucial feature of this form of the effective action is that it is local, in the sense 
that the form factors G n have an expansion in positive powers of the D%y This depends 
on making the choice of basis operators above, in contrast to the original form quoted in 



ref. 33 



The values of G n (0,0,0) for n = 0, . . .3 reproduce the Drummond-Hathrell results. 
For these, we have G n (0, 0, 0) = —\~ <?n(0, 0, 0), where: 



a = -~-sn(0,0,0) 

Z IT 



1 



a 



1 a 
2n 



03(0,0,0) 



144 7T 

1 a 
360 ^ 



b = 
d = 



la, 

02(0,0,0) 

I 7T 



1 a 

2 7T 



90 (0) 



13 a 

360tt 
1 a 

30 ?r 



(4.10) 



Explicit analytic forms for all the form factors G n are known and are given in detail in 
ref.fl8|]. Here, we simply quote the expressions for the form factors relevant for Ricci flat 



spaces, gz(xi,X2-,xz) and #9(3:1,3:2,3:3). Moreover, since terms involving derivatives acting 

A 2 

on the curvature are identified as higher order in j% in the present context, we restrict to 
the special = 0. Then, from ref.[P-8|, we find: 



fir 3 (0, £2, £3) 



and 



- F(0,x 2 ,x 3 )^ 4 + (l + ^(x 2 + x 3 ))(x 2 + x 3 ] 



+ /(* 2 )i 



3x 2 - x 3 1 



^3) 1- 



3x 



+ 



1/3:3 X 2 

-1 + 1 h x 2 + x 3 

2 \x 2 x 3 



2(x 2 + x 3 ) 2 
A 

2(x 2 + x 3 ) 2 
A 

1/1 1 

+ ^(— + — 

12 Va^ x 3 



1/1 1 



A\x 2 



,1 q 



(4.11) 



5r 9 (0,X 2 ,X 3 ) 



F(0, x 2 , x 3 ) — [4 + x 2 + x 3 ] 



+ f(x 2 ) 



A 2 {X2 



,1' q 



+ fM ~^ X 2 



Xc 



2x 2 



x 2 J 



2 



1 2 



(4.12) 
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where A = (x 2 — x 3 ) 2 . It can be checked that all the inverse powers of x 2 and x 3 cancel 

leaving a finite X2 = x 3 = limit, as required for the form factors to be local. Here, 

.1 



da e~ ail - a)3 



and 



F(xi,x 2 ,x 3 ) 



d 3 a 8(1 — a\ — a 2 — a 3 ) e 



■01-1012x3 — a>20t3X-L — ot3a-iX2 



a>0 



It will also be useful to note the simpler expression 

.1 



F(0,x 2 ,x 3 ) 



1 



X2 - X 3 J 



da — 

a 



t(l — a)x s 



a)x2 



(4.13) 
(4.14) 

(4.15) 



The behaviour of these functions is illustrated in the following plots of 173(0, x 2 , x 3 ) and 
(79(0, x 2 , x 3 ). Notice that along the diagonals x 2 = x 3 , both functions tend asymptotically 
to zero. However, if one argument is set to zero, then the functions may tend to a finite 
limit. The values at the origin are #3(0, 0, 0) = and (79(0, 0, 0) = 




Fig. 2 3D plots of 93(0, a?2, £3) over different ranges. 
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Fig. 3 Graphs of 93(0, X2, £3) along the X2 or X3 axes (left) and the diagonal xi = £3 (right). 
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Fig. 4 3D plots of 99(0, X2, £3) over different ranges. 




Fig. 5 Graphs of 39(0, X2, £3) along the X2 axis (top left), the diagonal X2 = £3 (top right) and the X3 
axis (lower). 
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5. Geometric Optics with the New Effective Action 



In this section, we apply the geometric optics methods introduced in section 3 to the 
effective action ( |4.7|) . In order to simplify the following discussion, we restrict here to the 
special case of Ricci flat spacetimes. This will be sufficient to extract the most important 
information, viz. the influence of the purely gravitational Weyl tensor contributions to high 
frequency photon propagation. 

It is easiest to present our results by building up from simplified cases@. So we consider 
first the modifications to the equation of motion = from the G3 and Gg terms 
neglecting their non-trivial derivative dependence. We then find: 



-4G 3 (0,0,0)1?/ap^ Ap + G9 (0,0,0) D 2 (R^ Xp D»F x ?) 



(5.1) 



Here, we have discarded all the 0(a) terms involving D^F^ for the reason explained 
in section 3, together with some terms involving derivatives of the curvature which are 
always suppressed. Notice that we are free to interchange covariant derivatives at will, 

A 2 

since a commutator produces a further power of curvature and therefore a further 0(j%) 
suppression. We have also used the Bianchi identities for the curvature and omitted terms 
involving the Ricci tensor. The important observation, however, is that the only effect on 
the equation of motion of the derivatives in the structure of the Gg term in the effective 
action is to produce an extra D 2 in the form factor. 

We now come to the implementation of the geometric optics approximation taking 
into account the non-trivial D 2 dependence in the form factors. Again, we illustrate this 
with a simple example. Consider a term 

l - [ dx^g R^x P F^D 2 F Xp (5.2) 



m 



in the effective action. This gives the following contribution to the equation of motion for 



D„F^: 



in 



-2(D 2 R^ Xp )D»F Xp - A(D a R^ Xp )D°D»F Xp - 4R^ Xp D 2 D»F Xp 



(5.3) 



Compared with the basic G3 term above, the relative orders of these terms are as follows. 

A 2 

The first is suppressed by O(-p-), since it involves extra derivatives of the curvature. The 
second involves an extra derivative of the curvature tensor but a compensating extra power 
of k from the associated derivative acting on the field strength, so overall this term is of 



5 Some early related analysis can be found in ref.[ 
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relative order O(jf^). In the low frequency region A > A c considered in M this would be 
neglected, but here we are interested in extending the results to A < A c . At first glance, 
the third term appears to be dominant because of the extra two powers of k coming from 
D 2 acting on F Xp . However, this appears contracted as k 2 which is not large but rather 
zero at leading order, and 0{a) in the full theory. So in fact this term must be discarded 
for consistency with the perturbative expansion in a. 

Returning to the second term, after substituting the geometric optics ansatz for F Ap , 
we find that we recover the same structural form for the light cone condition but with a 
factor now involving -^k ■ DR^ v \ p . It is clear that after we generalise from ( |5.2| ) to the 
complete form factors, these terms will sum up and we will find the light cone modified by 
functions of k ■ D acting on the curvature. 

A 2 

First, consider the order of magnitude of these corrections, O(jj^). This is interesting 
because it is closely related to the condition for direct observability of the Drummond- 
Hathrell effect. Consider differently polarised light propagating with a velocity difference 
of O (a -£§■), as predicted by eq.( fL2|) , over a time L characteristic of the spacetime curvature. 
This produces the biggest length difference between the rays which can be realised in the 
spacetime. To be observable, this should (as a rough order of magnitude) be greater than 
the wavelength A. We therefore arrive at the following criterion for direct observability: 

A 2 

a-%L > A (5.4) 

that is 

a AL > 1 (5 - 5) 
We see that to access the frequency range where superluminal effects could in principle 

A 2 

be observable, we need to satisfy the criterion jf- 3> 1. But since, as we have just shown, 
this is the parameter governing the corrections to the light cone from the new form factor 
terms in the effective action, we find that these terms cannot be neglected as perturbatively 
small but must instead be summed to all orders. 

Implementing this strategy, it can now be shown that the final light cone condition 
following from the equations of motion derived from ( f4.7| ) is: 



,2 2 

k 2 + 



m 2 



' / 2ik-D \ , „ / 2ik-D 
4G 3 0, =— ,0) -2ik-DGg[0, 5— ,0 



R^\pk^k x a v a p = (5.6) 



In the derivation, we have used the symmetry of G^{x\,X2,x^) under X2 <-> £3 (not present 

A 2 

in Gg of course), and omitted the O(j^) suppressed terms from the operator D 2 ^ acting 
on R pu \p. 
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Eq. (|5.6| ) is therefore the generalisation of the light cone (O), in the Ricci flat case 



where we extend the effective action to include the non-trivial form factor operators G n . 
In the light cone condition, these form factors reduce to single-variable functions of the 
operator k ■ D acting on the curvature tensor R pv \ p . Although the effective parameter 

1 A 2 

-^2 A; ■ D is O(jj^) and therefore not small in the region of interest, knowledge of the 
analytic expressions for the form factors should now enable these corrections to be exactly 
summed. We claim that this is the correct treatment of the equation of motion with extra 
derivatives in the perturbative corrections in a self-consistent geometric optics expansion. 

In the next section, we study the condition (|5.6|) in more detail and discuss what we 
can learn from it about the high frequency behaviour of the phase velocity v p h_- Before that, 
we make a few comments on the relation of our result to the analysis of Khriplovich |3T]] . 



In a very interesting contribution to the debate on dispersion and the Drummond-Hathrell 
mechanism, Khriplovich considered the general structure of the three-particle vertex for 
the scattering of an on-shell photon by a graviton. The 3 possible Lorentz structures are 
shown to be equivalent to the 3 terms in the Drummond-Hathrell action, acted on by form 
factors which are functions of the invariant momenta for the photon and graviton legs. 
The essential observation of Khriplovich is that for on-shell (i.e. k 2 = 0) photons, the form 
factors reduce to constants, given by the coefficients a, 6, c of section 3. The conclusion 
is then that there is no dispersion: the light cone shift at uj = is unchanged for all 
frequencies. 

Clearly, this differs from our conclusions here. The reason is not entirely clear. How- 



ever, the argument in |31| for the triviality of the form factors in the vertex looks remarkably 
similar to the discussion below ( |5.3|) . In our case, we also had form factors which were 
functions of D 2 acting on the electromagnetic field, and applying geometric optics these at 
first sight would appear to give contributions only of 0{k 2 ) which would vanish. However, 
the more complete analysis showed that they nevertheless can give rise to contributions 
involving k ■ D R which do not vanish. The most likely explanation for the discrepancy 



between our results and ref. [31] is probably the omission of the analogous terms in the 
analysis of the photon-graviton scattering vertices. 



6. High Frequency Photon Propagation 

The result of the previous section is the modified light cone formula for Ricci flat 
spacetimes: 

k 2 + g(^^) C^ Xp k»k x a»a p = (6.1) 
m \ m / 
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where G(x) is the known function: 

G(x) = -\- f°° - e~ is (is) g(sx) (6.2) 

2 7T Jq S 

where 

g(x) = 4(73(0, x, 0) + x g 9 (0,x, 0) (6.3) 

In terms of the functions f(x) and F(xi,x%, x^) defined in the last section, we can write 
the explicit form for g(x) as: 

16 ^ 11 5 1 

g(x) = -^F(0,x,0) + ^f(x) + {- + -)f'(x) + ^ + - (6.4) 
x z x A 2 x x z 15 

This is plotted below. Analytically, we can show g(0) = ensuring agreement with 
the Drummond-Hathrell coefficient c, while numerically we find that g(x) approaches an 
asymptotic value of 0.067 as x — > oo. 

.07 



0.06 
. 05 
.04 
0.03 



10 20 30 40 50 

Fig. 6 Plot of the function g(x) which enters the formula for the modified light cone. 

With this explicit knowledge of the form factor, in principle we have control over the 
high frequency limit of the light cone and the phase velocity. However, the next difficulty 
is that in eq. fl6.1T ), G is a function of the operator k ■ D acting on the Weyl tensor. If 
C^vXp is an eigenfunction then G will reduce to a simple function of the eigenvalue and 
then Fig. 6 determines its asymptotic behaviour. However, in general this not true and 
consequently the interpretation of ( |6.1| ) is far from obvious. 

At this point, the problem is reduced to differential geometry. The encouraging feature 
of (|6.1|) is that the operator k ■ D simply describes the variation along a null geodesic. 
(Notice that because the second term is already 0(a), we can freely use the usual Maxwell 
relations for the quantities occurring there, e.g. the results k ■ Dk u = and k ■ Da u = 
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derived in section 3.) The question then becomes what is known about the derivative of 
the Weyl tensor along a null geodesic. 

As in our previous work, it is convenient to re-cast the light cone condition in Newman- 
Penrose formalism (see, e.g. refs.[^5j,[^6l for reviews). This involves introducing a null 



tetrad e a M based on a set of complex null vectors n^, m M , m^ 4 ) . The components of 
the Weyl tensor in this basis are denoted by five complex scalars ^0, . . . , ^4. For example, 
* = -C abcd £ a m b £ c m d and # 4 = -C abcd n a rh b n c rh d . 

Consider propagation along the null direction i^, i.e. choose /c M = . The two 
spacelike polarisation vectors a M and are related to the m M and fh^ null vectors by 
m M = A= (a M + m 1 * = -j= [a^ — In this case, the light cone condition can be 

simply written as 

k 2 + 4 G ft^)(*o + <) = (6.5) 

We have searched the relativity literature for theorems on the behaviour of I^D^q 
without finding any results of general validity, although it seems plausible to us that some 
general properties may exist. (Notice, for example, that naturally enough this is not one of 
the combinations that are constrained by the Bianchi identities in Newman-Penrose form, 
as displayed for example in ref.|3(|, ch.l, eq.(321).) To try to build some intuition, we 
have therefore looked at particular cases. The most interesting is the example of photon 



propagation in the Bondi-Sachs metric [ p7| , |38| which we recently studied in detail 



The Bondi-Sachs metric describes the gravitational radiation from an isolated source. 
The metric is 

ds 2 = -Wdu 2 - 2e 2l3 dudr + r 2 h l3 (dx % - U l du)(dx j - U j du) (6.6) 

where 

h lJ dx l dx J = \{e 21 + e 25 )d6 2 + 2 sinh( 7 - 5) sin6W0# + \{e~ 21 + e~ 2S ) sin 2 9d(f> 2 (6.7) 

The metric is valid in the vicinity of future null infinity X + . The family of hypersurfaces 
u = const are null, i.e. g^d^ud^u = 0. Their normal vector satisfies 

e„ = d^u i 2 = 0, e^D^r = o (6.8) 

The curves with tangent vector £ M are therefore null geodesies; the coordinate r is a radial 
parameter along these rays and is identified as the luminosity distance. 

The six independent functions W. /3, 7, 5, U l characterising the metric have expansions 
in - in the asymptotic region near X + , the coefficients of which describe the various features 
of the gravitational radiation. (See [|J for a brief summary.) 
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In the low frequency limit, the light cone is given by the simple formula (|1.2j ). The 
velocity shift is quite different for the case of outgoing and incoming photons [|J. For 
outgoing photons, /c M = £ M , and the light cone is 



while for incoming photons, = n M , 

Now, it is a special feature of the Bondi-Sachs spacetime that the absolute derivatives 
of each of the Weyl scalars V&Oj • • • » ^4 along the ray direction vanishes, i.e. ■ • ■ ■> ^4 



are parallel transported along the rays p8| , p5| . In this case, therefore, we have in particular: 



£-D ^ = £-D^ 4 = (6.11) 

but there is no equivalent simple result for either n ■ D f 4 or n ■ D f o- 

These results can be easily checked for the simpler related example of a weak-field 
plane gravitational wave [|J. In this case, the first identity is trivial since \I/o = 0, but 
in particular we can confirm that n ■ D \l/4 7^ and ^4 is not an eigenfunction. The 
important result i ■ D ^ = therefore appears to be a very special property of the 
Bondi-Sachs metric and not an example of a general theorem on derivatives of the Weyl 
tensor. 

Although it seems to be a special case, ( |6.11|) is nevertheless important and leads to 
a remarkable conclusion. The full light cone condition ( |6.5| ) applied to outgoing photons 
in the Bondi-Sachs spacetime now reduces to 

k 2 + ^G(0)(* + %) = (6.12) 

since £ ■ D^o = 0. In other words, the low frequency Drummond-Hathrell prediction of a 
superluminal phase velocity t> p h(0) is exact for all frequencies. There is no dispersion, and 
the wavefront velocity u p h(oo) is indeed superluminal. 

This is potentially a very strong result. According to the analysis presented here, we 
have found at least one example in which the wavefront truly propagates with superluminal 
velocity. Quantum effects have indeed shifted the light cone into the geometric spacelike 
region, with all the implications that brings for causality. 
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7. Propagation in Background Magnetic Fields 



Before accepting the results of the last section as definitive, however, it is instructive 
to compare what we have done in this paper for the case of a background gravitational 
field to previously known results on the propagation of light in a background magnetic 
field. 

The refractive index for photons moving transverse to a homogeneous magnetic field 
B has been calculated explicitly as a function of frequency in two papers by Tsai and Erber 
3^,40]. They derive an effective action 



dx 



(7.1) 



where M^„ is a differential operator acting on the electromagnetic field A v . Writing the 
corresponding equation of motion and making the geometric optics ansatz described in 
section 3, we find the light cone condition 

k 2 - M^(k)a»a v = (7.2) 

In terms of the refractive index n(uj) = c— = ct> p h(w) -1 , this implies 

n(u) = 1 - -^M^(k)aW (7.3) 

Denoting M^ v a^a v by My, Mj_ for the two polarisations, the complete expression for 
the birefringent refractive index iM: 

n\\,±(u) = 1 " ^2 M IU 

a (eB\ 2 f 1 f°° -is(i+s 2 n 2 p(u,z)) 

(7.4) 

where z = and O = — . The functions iV and P are given by: 

x 1 ( cos zu — cos z l—u 2 \ 1 , A 2 \ 1\ /B r > 

P(u,z) = -^(—-—. — = — (1-u 2 ) + 0(z 2 ) 7.5 

z z V Zzsmz 4 / 12 



6 The extra power of (is) in the prefactor relative to eq.( |6.5| ) arises because the lowest order 
terms in the effective action contributing to ( |7.1|) are of 4th order in the background fields, i.e. 
0(F 4 ), compared to the 3rd order terms of O(RFF) in the gravitational case. For background 
magnetic fields, the analogue of the Drummond-Hathrell low frequency action is the familiar 
Euler-Heisenberg action (see, e.g. ref.0]). 
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and 

Nn 



cot z 
2z 



o usmzu\ cos zu 
1 - u 2 + — + 



sin z 



2zsinz 



N 



± 



z cos zu zu cot z sin zu zicos zu — cos z) 



2 sin z 



2 sinz 



sin 2 z 



i(i-« 2 )(i + ^ 2 )+o(^ 2 ; 



In the weak field, low frequency limit, we can disregard the function P and consider 
only the lowest term in the expansion of iV in powers of z. This reproduces the well-known 
results (see also [ffljpl): 



n \\,± 



a feB\ 2 
,o 1 + — — 
47T V m / 



14 8 
L45|i' 45±J 



(7.7) 



The weak field, high frequency limit is analysed in ref.[^J. It is shown that 

a feB\ 2 r -i . , 

n hl _ ~^oo 1-— ^— j [ C ||,cJO 3 (7.8) 

where the numerical constants are [c||,cj_] = 3 1^ v /7rr(|) 2 (r(|)) 1 [3||,2j_]. 



The complete function n(uS) is sketched in Fig. 7. It shows exactly the features found 
in the simple absorption model described in section 2. 



n(w) 




Fig. 7 Sketch of the frequency dependence of the refractive index n(uj) for light propagating in a background 
magnetic field. The crossover point is at £1 ~ 1. 

In particular, the phase velocity v p h(uj) begins less than 1 at low frequencies, showing 
birefringence but conventional subluminal behaviour. In the high frequency limit, however, 

4 

the phase velocity approaches c from the superluminal side with a &~z behaviour. 
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All this is by now standard and in line with our expectations. What we are interested 
in here, however, is the comparison between this calculation of v p h(u)) via the effective 
action ( |7.1| ) and the gravitational calculation of sections 5 and 6, in particular the origin 
of the non-analytic high frequency behaviour ( |7.8| ). 

The important observation is that the exponential term involving 2 P is crucial in 
obtaining the high frequency limit. However, if we made a literal expansion of the effective 
action in powers of the magnetic field (eS/m 2 ), this term would be regarded as higher 
order and discarded. It is important however, because it involves the product of the field 
and the frequency, (eBui/m 3 ), and this is not small in the interesting high frequency region. 
The exponent in the effective action is vital in describing the high frequency propagation. 

Now compare with our construction of the effective action in a background gravita- 
tional field. The action ( [4.71) is obtained by expanding in the curvature and keeping only 
terms of O(R). This appears to be analogous to the expansion of iV in (|7.4j ). We have 
found a non-trivial extension of the zero- frequency Drummond-Hathrell result by taking 
into account derivatives of the curvature. These are essential in the gravitational case 
because there is no change in the light cone for a constant curvature metric (since this 
is totally isotropic), whereas there is an interesting birefringent effect even for a constant 
magnetic field. By restricting the effective action rigorously to terms of first order in the 
curvature, however, we seem to have missed the analogue of the exponent terms in ( |7.4| ) , 
which would be characterised by the not necessarily small parameter (Ru/m 3 ). 

The conclusion of this comparison is therefore discouraging. It seems that despite its 
complexity, the effective action (|4.7| ) may still not be sufficiently general to encode the 
high frequency behaviour of photon propagation. 



8. Conclusions 

In this final section, we attempt to synthesise what has been learnt from this investiga- 
tion and identify what further work is required for a complete resolution of the dispersion 
problem for the propagation of photons in gravitational fields. 

In view of our experience with the background magnetic field problem, it is natural to 
assume that the full effective action governing propagation in a background gravitational 
field takes an analogous form and the modified light cone can be written heuristically as 

fc2 + a f°° da e -is(i+sWv(s,R)) = (81) 

Wo s 

where both Af and V can be expanded in powers of curvature, and derivatives of curvature, 
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with appropriate powers of s. The frequency dependent factor would be Q ~ -^--^ ~ 

A 3 

0(^3-), where 'it!' denotes some generic curvature component. If this is true, then an 
expansion of the effective action to first order in O(^) would not be sensitive to the V 
term in the exponent. The Drummond-Hathrell action would correspond to the leading 
order term in the expansion of J\f(s, R) in powers of -M? neglecting derivatives, while our 
improved effective action sums up all orders in derivatives while retaining the restriction 
to leading order in curvature. 

The omission of the V term would be justified only in the limit of small O, i.e. for 

A 3 

<C 1. Neglecting this therefore prevents us from accessing the genuinely high frequency 
limit A — > needed to find the asymptotic limit i> p h(oo) of the phase velocity. If eq. (|SJ]) is 
indeed on the right lines, it also looks inevitable that for high frequencies (large O) the 
rapid variation in the exponent will drive the entire heat kernel integral to zero, restoring 
the usual light cone k 2 = as uj — > 00. However, this is simply to assert that in this respect 
the gravitational problem behaves in the same way as the magnetic field background, and 
while this seems plausible in relation to the ui — > 00 limit we should be cautious: systematic 
cancellations or special identities (c.f. section 6) could change the picture, and we should 
remember that it was not foreseen that the low frequency behaviour of f p h(0) would differ 
so radically from other background field calculations as to produce superluminal velocities. 
So although the available evidence seems to point strongly towards a restoration of the 
usual light cone in the high frequency limit, this inference should be made with some 
caution. 

Unfortunately, the quantum field theoretic calculation required to settle the issue by 
evaluating the l V type correction to the exponent in ( |8.X[ ) looks difficult in general, at least 
comparable to the evaluation of the effective action 'A/ - ' terms in section 4. However, if we 
are only interested in the to —>■ 00 limit, it may be sufficient just to perform some leading 
order type of calculation to establish the essential (l + s 2 Q 2 V) structure of the correction. 
Another possible simplification, which would be interesting in its own right, would be 
to reformulate the effective action calculation from the outset in Newman-Penrose form, 
thereby removing the plethora of indices arising in calculations involving the curvature 
tensors. Simple cases, such as black holes, in which only one Weyl scalar is non-vanishing 
might prove particularly tractable. 

However, even if this picture is correct and the light cone is eventually driven back to 
k 2 = in the high frequency limit, the analysis in this paper still represents an important 
extension of the domain of validity of the superluminal velocity prediction of Drummond 
and Hathrell. Recall from section 5 that the constraint on the frequency for which the 

A 2 

superluminal effect is in principle observable is yfj ^> 1. Obviously this was not satisfied 
by the original uj ~ derivation. However, our extension based on the 'all orders in 
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derivatives effective action does satisfy this constraint. Combining with the restriction 

A 3 

Y^r <C 1 in which the neglect of the V type corrections is justified, we see that there is a 
frequency range 

T » i» % ^ 



for which our expression fl6.1| ) 



k 2 + 4r Gi 2 ^-) C^ Xp k»k x a»a p = (8.3) 



for the modified light cone is valid and predicts observable effects. 

Since this formula allows superluminal corrections to the light cone, we conclude that 
superluminal propagation has indeed been established as an observable phenomenon even 
if, as seems likely, causality turns out to respected through the restoration of the standard 
light cone k 2 = in the asymptotic high frequency limit. 
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